GLOBAL EXISTENCE AND UNIQUENESS OF WEAK 
SOLUTIONS OF 3-D EULER EQUATIONS WITH HELICAL 
SYMMETRY IN THE ABSENCE OF VORTICITY STRETCHING 



BORIS ETTINGER AND EDRISS S. TITI 



Abstract. We prove uniqueness and existence of the weak solutions of Euler 
equations with helical symmetry, with initial vorticity in L°° under "no vortic- 
ity stretching" geometric constraint. Our article follows the argument of the 
seminal work of Yudovich. We adjust the argument to resolve the difficulties 
which are specific to the helical symmetry. 



1. Introduction 

Ideal incompressible homogeneous fluid of density po a-nd confined in three- 
dimensional domain D C K-^ is governed by the Euler equations: 

(1.1a) ^ + (u. V)u = --Vp + F, 

dt Pq 

(1.1b) V-u = 0, 

supplemented with the no-normal flow boundary conditions 

(1.2) u • n = 0, on 9D, where n is the normal vector to 9D, 

and initial velocity uo(x). u : D x [0, T) is the velocity field, p : D x [0, T) ^ M 

is the pressure, determined by the incompressibility condition and F : D x [0, T) ^ 
M.^ is the given external body forcing term. We will consider constant density 
Po = 1- 

In this article, we will investigate the solutions of equations (1.1), which are 
invariant under a helical symmetry group C^. The group C is a one-parameter 
group of isometries of M.^ 

(1.3) = {S-p : ^ M^Ip e M}. 

The transformation Sp {S stands for "screw motion") is defined by: 

(x\ fx cos p + y sin p 
y\ = -xsinp-l-ycosp 
zj \ Z + Kp 

where k is a fixed nonzero constant length scale. In fact, Sp is a superposition of 
a simultaneous rotation around the z-axis with a translation along the z-axis. The 
symmetry lines (orbits of C^) are concentric helices. We call the solutions, and 
more generally functions, which are invariant under - "helical" . Since the Euler 
equations in M.^ are invariant under isometries, then (under mild assumptions of 
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uniqueness) solving the Euler equations in a domain which is invariant under heUcal 
symmetry with a heHcal initial condition and a helical body forcing will give rise 
to a solution, which is helical for the whole interval of time of it's existence. 

Observe that S2tt is a translation by 27rK in the z direction. Therefore, helical 
symmetry imposes a periodic boundary conditions in the z direction. We will also 
assume that the physical domain D C is bounded in the x and y directions, thus 
imposing a no-normal flow. We say that D is bounded in an (infinite) cylinder along 
the z-axis, which has a finite radius. Since D is invariant under G"^, it's boundary 
5ID) can be thought as being the union IJ SpC, where C is a closed planar curve. 

The difficulty to establish the global regularity for the solutions of the 3D Euler 
equations can be best appreciated, when one examines the evolution of the vorticity 
field = V A u, where VA is the curl (rotor) operator. Taking the curl of both 
sides of equation (1.1a) we get: 

(1.5) ^ + (u- V)n + (n- V)u = VAF. 

The last term of the left-hand side, ($7 • V)u, is called the vorticity stretching 
term. This term is the main obstacle to achicivc global in time regularity of the 
three-dimensional Euler Equations, (sec [3], [6], [17] for the latest discussions of this 
question). The difficulty remains after imposing the helical symmetry on the solu- 
tion, since helical flows can undergo nontrivial vorticity stretching. 

We will therefore include an additional requirement. We will demand that the 
velocity field u = {ux,Uy, Uz)'^ , where Mx, Uy, Uz are components of the vector field 
in the basic directions, obeys the following constraint: 

(1.6) yux — xuy + KUz = 0. 

This condition is an orthogonality of the velocity field to the symmetry lines of the 
group C^. This condition together with the assumption of the helical symmetry 
lead to vanishing of the vorticity stretching term. We will prove in Section 2 that 
under these conditions, the vorticity ft is directed along the symmetry lines and 
it's magnitude, up to a normalization is transported by the flow. This control of 
the L°° norm of the vorticity is the key to our argument of global existence and 
uniqueness. It is consistent with the celebrated result of Beal-Kato-Majda [4], and 
the works of Yudovich [23] and Ukhovskii and Yudovich [19] for the 2D case and 
axi-symmetric flow (without swirl), respectively. 

Our article is inspired by the seminal work of Yudovich [23] , who proved existence 
and uniqueness for a certain class of weak solutions in two space dimensions. We 
adapt his ideas: the stream function weak formulation of the problem and the 
elliptic regularity to the case of helical flows. We follow closely his article, especially 
in Section 4, where we prove uniqueness. A different route would have been to 
extend the ideas of Bardos [2] , adding viscosity with artiflcial boundary conditions 
and passing to the limit of viscosity going to zero. In such case, the equation 
with the added viscosity is not the Navier-Stokes equation, because of the different 
boundary conditions and the limit also has a non-physical boundary conditions. 

In his work, Yudovich [23] proved global existence and uniqueness of the solutions 
of the two-dimensional Euler equations, whose initial vorticity belongs to the space 
of essentially bounded functions, L°°. Uniqueness of the solutions was extended 
for a wider classes of functions in the work Yudovich himself [24] to the class 
of functions, which are not bounded but whose norms grow "slowly enough". 
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Uniqueness and global existence was also proven by Vishik [20] for vorticity in 
the Besov-like spaces. Global existence of solutions of the two-dimensional Euler 
equations was proven for initial vorticity CIq G by Majda and DiPerna [9] and 
also for solutions, whose initial vorticity is a positive Radon measure by Delort [8] . 
No uniqueness is known in these cases. 

The helical flows fall within a class of "two-and-a-half dimensional flows, namely 
flows in a three-dimensional domains with a certain continuous spatial symmetry. 
The most heavily investigated within this class are the axi-symmetric flows, which 
are invariant under a rotation around a certain axis of symmetry. The viscous 
axi-symmetric flows^ were analyzed by Ladyzhenskaya [15], where the author had 
to assume that the domain of the flow does not contain a cylinder of a positive 
radius around the axis of symmetry, in order to prove global well-posedness. The 
inviscid axi-symmetric flows were the subject of Ukhovski and Yudovich [19], the 
authors had to assume that the vorticity vanishes rapidly enough near the axis of 
symmetry, namely j e L°° to conclude global existence and uniqueness of the weak 
solutions. Both of these articles assume that the azimuthal, "swirl" component ue 
of the velocity is zero. 

The assumption of zero azimuthal components for the axi-symmetric flows means 
that the vorticity stretching term in (1.5) is zero. There were several attempts to 
analyze the axi-symmetric flows with non-trivial vorticity stretching term. Chae 
and Immanuilov [5] proved that solutions exist globally in time and are unique 
for a "generic" class of solutions. Moreover, Hou and Li [13] constructed recently 
a family of infinite energy solutions that can undergo a very rapid and intense 
vorticity stretching yet do not develop a singularity and exist globally in time. 

The global in time existence and uniqueness of viscous helical flows for all helical 
initial data was established by Malialov et al [16] without any restrictions, such 
as zero vorticity stretching, on the velocity field. Inviscid helical flows were the 
subject of Dutrifoy [10]. The author of [10] proves global existence and uniqueness 
of the classical solutions whose velocity field is orthogonal to the symmetry lines, 
by essentially using the fact that the vorticity is transported by such a flow and 
applying the Theorem of Real, Kato and Majda [4]. In this paper, we extend the 
conclusions of [10] to a weaker class of solutions. 

The helical symmetry introduces a complication, as opposed to planar and axi- 
symmetrical cases, because one cannot flnd an orthogonal coordinates, such that 
one of them runs along the symmetry lines. Nevertheless, we were able to reduce 
the Euler equations in our problem to effective equations for two-dimensional ve- 
locity on a two-dimensional domain. Using this reduction, we produce an efl'ective 
vorticity-stream function formulation, which allows us to proceed along the lines of 
Yudovich's work [23]. A similar reduction was proposed in [10, Section 3], but was 
not investigated to the full extent of it's consequences. 

This article is organized as follows. In Sec;tion 2, we lay down the geometric 
framework and obtain the consequences of our assumptions to classical solutions 
of Euler Equations, with this we perform the reduction to the two-dimensional 
formulation. In Section 3 we introduce the appropriate functional-analytic setting 
and formulation of the weak problem. In Sections 4 and 5, we prove the uniqueness 
and the existence of the solutions of the weak problem, respectively. 



In this article, we will speaJj always about incompressible flows. 
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2. Geometrical Setting and Stream Function Formulation 

In this section, we investigate the geometry of helical symmetry and of the or- 
thogonality assumption, Equation 1.6. Wc; then rewrite the problem of the heli- 
cal velocity field, which obeys the orthogonality assumption as a problem of two- 
dimensional velocity field on a two-dimensional domain which obeys a set of effective 
equations. We introduce the effective vorticity stream function formulation, which 
will be the basis of the analysis further in the article. 

A note on our notation. We will usually specify vector and a vector field in a 
boldface e.g. u and it's coordinates in a normal font, with indices x, y, z or 1, 2, 3 
e.g. u = {ux,Uy,Uz)^ or u = {ui,U2,U3)-'^ . We denote by (•,•) the Euclidean inner 
product in M^. It is also used to denote the pointwise inner product of two vector 
fields. 

2.1. Invariance under helical symmetry. Denote by Rp the rotation transfor- 
mation around the z axis by an angle of p radians: 

(2.1) Rp = 

With this notation, the action of on K"* becomes: 

(0 \ / X cos p + y sin p 
= -xsinp + ycosp 
Kp) \ Z + Kp 

We solve our equation in a domain B. Since we look for solutions, which obey 

boundary conditions imposed by D, and since these solutions have helical symmetry, 
then D must also be invariant under a helical symmetry, namely be a union of 
helices. 

Definition 2.1 (Helical Domain). A domain Hi C is helical, if it is invariant 
under the action ofC^, namely 

(2.3) SpB = P, Vp e R. 

For the rest of this article, we will assume the following assumption about the 
domain D: 

• D is helical, 

• ID) is simply connected. 

• D is contained in a cylinder of finite radius around the z axis. 

• B has a twice diff'erentiable boundary. 

The sole purpose of the simply connectedness assumption is to simplify our proof, 

it can be removed to yield a more general result. 

The action of G"^ on a scalar function is by composition: 

(2.4) /^''(x)=/(5_px). 

The action of on vector fields has to take care of the underlying rotation: 

(2.5) v^''(x) =i?pv(5_px). 

The following vector field will have an important role throughout this article: 

(2.6) ^ = yx — xy + kz. 
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The vector field ^ is tlie field of tangents of the symmetry lines of G". 

Definition 2.2 (Helical Function). A (scalar) Junction f is called helical, 

if 

(2.7) /(5,x)=/(x), VpeM. 

Namely, it is invariant under the action of . 

Claim 2.3. A continuously differentiable function f -.M.^ is helical if and only 

if 

(2-8) ^i-^|+'^|f = (^>^)=0- 

Proof. Assuming that the function / is helical, differentiate the relation 

(2.9) /(5px) = /(x) 

with respect to p. The right-hand side is zero, the left-hand side is 

[f{x cos p + y sin p, —x sinp + y cos p,z + np)] 

^^■"^^^ df df df 

= {—xsinp + y cos p)— ^ — h {—x cosp — j/sinp)— ^ — h k-;^-. 

ox oy oz 

Now compute at p = to get the desired relation. 

To prove the converse of the statement, let x e R^, denote c = (ci, 02,03)^ : 
R ^ R3 by c(p) = SpK. Define a function 5 : M M by: 

(2.11) 9{p) = f{c{p)). 

By virtue of (2.8), the function g{p) satisfies the following differential equation: 

(2.12) g = c2(p)|^(c(p))-ci(p)|^(c(p))+Kg(c(p)) = 0, with ff(0) = /(x). 
Therefore, 

(2.13) g{p) = g{0), and consequently, /(5px) = /(x). 

□ 

We will also use the following notation: 

With this notation the last claim states that a differentiable function / is helical if 
and only if || = 0. 

Definition 2.4 (Helical vector field). The vector field v : 1— > is helical, if 
(2.15) v(6'px) = i?pv(x), Vp e M. 

Claim 2.5. A continuously differentiable vector field v = {vx,Vy,VzY' : — > R^ 
is helical, if and only if it obeys the following relations: 
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(2.16c) ^ = 0. 

Proof. Assume the relation (2.15) and differentiate it with respect to p at p = 0. 
In the right-hand side, only Rp depends on p and we have: 



p=0 




In the left-hand side, by differentiating we obtain a componentwise differentiation 
along ^. Comparing the two expression, we get the required equality. To prove the 
converse, we apply the existence and uniqueness theorem for a system of ordinary 
differential equations in a similar manner to the proof in the case of the scalar 
helical function in Claim 2.3. □ 



We will introduce the following notation 

■ 

p=0 



(2.18) n - '^^^ 




dp 

With this notation, the previous claim can be stated in the following manner: 
a continuously differentiable vector field v is helical, if and only if v obeys the 
following relation 

(2.19) 5 = ^^- 

Let us write the relation that helical vector fields obey, i.e. (2.16) in a more 
explicit manner: 

(2.20a) - ^ a7 + " 9^ = 

dVy dVy dVy 

(2-20b) y^--^+-^ = -^- 

Definition 2.6. Let \i,p he a solution of Euler equations then we call u,p a helical 

solution, if u is a helical vector field and p is a helical scalar function. 

2.2. The orthogonality condition. We will introduce the following notation, for 
a vector field v = {vx, Vy,VzY' , we will define a function by: 

(2.21) V^= [v.Cj = - XVy + KV^. 

Lemma 2.7. Let \i,p he a smooth helical solution of the Euler equations (1.1), 
then: 

Due 

(2.22) 

where is the material time derivative, defined by ^ = ^+Ux§§+Uy^+Uz^. 
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Proof. The full form of the three-dimensional Euler equations is: 
(2.23a) — + + uy— + = -- + Fx, 



(2.23b) 


dUy 


dUy 

OX 


duy 


oz 


dp 
dy 




(2.23c) 


at 


duz 
ox 




duz 


dp 
dz 





We take x times equation (2.23a) —y times equation (2.23b) and k times equation 
(2.23c). The right-hand side of (2.23) gives 



(2.24) 



d^ ^ ^ 



where we used the fact that p is a helical function. Prom the left-hand side of (2.23) 
we get: 



(2.25) 



du^ 

'df 


+ y 


dux 
dx 


dux 


dux' 
dz 




+ K 


duz 
dx 


duz 


duz 
dz 



— X 



dUy dUy dUy 

dx dy dz 



du^ 



dur 



du^ du^ 

_+Ux—+ UxUy + Uy^- UxUy + n,— 



By equating the two sides we conclude the lemma. 



□ 



Corrolary 2.8. Suppose uq, F are helical vector fields that give rise to a helical 
smooth solution u,p of 3D Euler equations. If F^= and Uq^ = then u^t) = 

oyt. 

Definition 2.9. We will say that a vector field v is orthogonal to the helices, if 
(2.26) = yvx - xvy + kv^ = 0. 

Definition 2.10 (Orthogonality Assumption). Throughout this article we will as- 
sume that the forcing field F and the velocity field of the solution of the Euler 
equations u are orthogonal to the helices in the sense of Definition 2.9. Namely, 



(2.27) 



Fr=0, 



(2.28) U|- = yux — xUy + kUz = 0, 

for interval of time for which the solution exists. 

Remark. By Lemma 2.7, it is enough to assume that F^ = for all times and 
|- = {£,, Uq) = in order to satisfy the previous definition. 

The following two lemmas are the direct consequences of the Orthogonality As- 
sumption (Definition 2.10) about the velocity field on it's vorticity J2 = V A u. 
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Lemma 2.11. Let u be a helical vector field, such that — 0. Denote ft 
V A u = {Qx, ^y, ^z)'^ , the vorticity of u then 

LU 



(2.29) 



n 



= (yOJ, —XLU, Kiu) 



for a helical scalar function w = fl^ = 



dUy 

dx 



dux 



Proof. Observe that if representation (2.29) holds true then since fi has zero diver- 
gence, we have 



du) du 



(2.30) y X—+K — 

ox oy oz 



du> dfl-r dfl^i dflz 



dx 



+ 



dy 



+ 



dz 



V • n = 0. 



Therefore by Claim 2.3, w will be a helical function. We compute the components 
of ft: Observe that equation (2.26) can be rewritten as: 



(2.31) 



1/ 

-y-yux + xuy). 

K 



Also, by using equations (2.20a) and (2.20b), we can conclude the following equal- 
ities: 



(2.32) 
(2.33) 



dux _ dux 
dz K dx 



dux _ , 

^ Qy Uy), 



dUy 1 , dUy dUy , 

dz K dx dy 



Therefore, we use (2.31) and (2.32) to obtain: 
(2.34) 



du^ 
dy 
1 " 

K 



dUy 19. ,1 

= -T- [-yux + xuy] — 

dz K dy K 

dux duy duy 

- - +Ux + y-^ 

dx 



dUy 
dx 



dUy 

dy 



dy dy 
and we use (2.31) and (2.33) to obtain: 
(2.35) 



- y. 


dUy 


dux 


K 


dx 


dy 



Oil — 



dux duz 
dz dx 



Uy - y- 



dux 
dx 



dux 
' dy 



dy 



Id, 



dux dux dux duy 
dx dy dx dx 



—x 


dUy 


dux 


K 


dx 


dy 



(2.36) 



di 



dux 



dx dy 



We see from the calculations above that 



(2.37) 



n = ^- 



dUy 

dx 



dux 
dy 



proving that VI has the representation (2.29) with uj = — = Clz- 



□ 
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Remark. Observe that the previous lemma also proves that ifF is a body forcing 
term, which obeys orthogonality condition, then the vorticity forcing term V A F 
obeys 

(2.38) VAF = ^(VAF);,. 

Lemma 2.12. The orthogonality condition has also direct implications to the dy- 
namics of the vorticity. Let Jl = V A u be the vorticity field of a smooth helical 
solution of the Euler equations which obey the Orthogonality Assumption (Defini- 
tion 2.10), then n obeys the following equation 

DU 1 

(2.39) — + -^7^7^u = VAF. 

Ut K 

The matrixTZ is defined in equation (2.18). 
Proof. From (1.5), the vorticity fl satisfies 

(2.40) ^ + ■ '^^^ + (n • V)u = V A F. 
By Lemma 2.11 , we have 

(2.41) ft^ = ^-flz. 

(2.42) fly = --Oz. 



K 



Therefore, 



{fl ■ V U = fl^— + fly— + fl,— 

OX oy oz 

^^■^^^ 1 / 5u ^ du ^ du\ 1^ du 1^ ^ 

= - y^z-^ a;S2^— + kS]^— = -flz—=i = -O^K-u, 

K \ ox oy oz J K K 

where we used the fact that u is a helical vector field and Claim 2.5. □ 

Corrolciry 2.13. The helical scalar function uj = flz, where fl = V A u, the 
vorticity field of a smooth solution of the Euler equations, obeys the equation 

(2.44) D.^djy_m_ 
^ ' Dt dx dy 

Remark. Observe from (2.29) that 

(2.45) ui = flz = \fl\ 



Therefore, the equivalent of the previous two lemmas in the axi-symmetric setting 
would be the fact that if the azimuthal component of the velocity is zero - uo = 

D 
Dt 



then ft = (JeO and (^) = (V A F)^. 



2.3. Reduction to the two-dimensional equations. The purpose of this sub- 
section is to setup the reduction of the three-dimensional problem to the problem 
with two components of velocity fields on a two dimensional domain. To achieve 
this goal, we rewrite the orthogonality condition (equation (2.26)) in the following 
form: 

(2.46) Uz = -(-ywx + xuy). 
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In the equations (2.20), which characterize the hehcal vector field, we move differ- 
entiation by z to the left-hand side: 

(2.47a) «_ = ^^_y_+^_ 



duy duy d 



X- 



(2.47b) ^^ = --^-ydx -dy 

(2.47c) .^=-.^+.^. 

oz ox oy 

Prom equations (2.46) and (2.47) it is clear that we can eliminate the Uz component 
from the velocity field and also all the differentiation by z. It means that we can 
rewrite our problem, as a problem for a two-dimensional vector field Uj-ed = {ux, Uy), 
with the fluid dynamics occuring in a two-dimensional domain of a form = 
D n {z = po}. We will make an arbitrary choice of po = 0. We define a subset of 
R2 by 

(2.48) Iio = {{x,yf\{x,y,Of€B}. 

We will solve the new equations on Dq and then we can recover the full vector field 
on D by the following algorithm: 

Input: 

(2.49) u^{x, y, 0; t),Uy{x, y, 0; t) on ©o x {0}. 
Step One: Recover the Uz component on Dq x {0} by: 

(2.50) Uz{x, y, 0; t) = -{-yux{x, y, 0; t) + xuy{x, y, 0; t)). 

K 

Step Two: Recover the velocity field on P from the velocity field on Pq by: 

(2.51) u{x,y,z;t) = R^u{S-±x;t), 
where x = {x, y, z)'^. 

Observe that in a similar manner, we eliminate differentiation by z from any helical 
function / by using equation (2.8): 

2.52 ^J=-yJ+x^. 

oz ox oy 

Definition 2.14. For a function (j) : Dq ^ R denote by 4> the helical function, 
which extends (j) to D; 

(2.53) (i){x,y,z) ^ (l){S-i.{x,y,z)^). 

For a helical function ij) •.'O ^'U. denote it's contraction to Dq by 'ip. 



(2.54) 



V'(a;,t/) = V(a;,y,0). 
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2.4. Effective vorticity-stream function formulation for helical solutions 
of the Euler Equations. We will now introduce the effective stream function 
formulation on the domain Dq. We rewrite the incompressibility condition: 
dUx duy ^ du;^ _^ duy y du^ _^ x du^ 

dx dy dz dx dy k dx k dy 
duy yd X d 

^^■^^^ = 97 + a7 " ^"^"^ + '^"''^ + + ^"-^ 

= + y^)"^^ - ^y^vi + ^'^^^'^^ +x'^)uy- xyu^]. 

We used equation (2.20c) in the first equality and then equation (2.46) in the second 
equality. Therefore the incompressibility for a reduced velocity field is: 

(2.56) ^^[^'^^ ^^^"^ ~ ^y^y^ /t2 ^[^'^^ ^^^"^ ~ ^^"^^ ^ ^ 

We introduce the function i/' : Dq — > R by 
(2.57a) ^ = -1^ [{k^ + y'^)ux - xyuy] , 

(2.57b) [-xyux + {i^^ + x'^)uy] . 

Let us rewrite the relation between tp and Ux,Uy: 

(2.58) fiV^f 2 + 

Now, we can invert the previous equation to obtain: 
(2 59) ( ""A = ^ ( " 

^ ' \UyJ k'^ + y2 + y2 

Remark. In fact, the function tp has an intrinsic meaning to the full three-dimensional 
problem. Consider the three-dimensional stream function of the problem, defined by 
V A = u. Unlike the planar and axi-symmetric problem, it is no longer true that 
* = S.ip, where tp is the helical extension of ip in the sense of Definition 2.1J^. But 
ip can be recovered as a component of ^ in a certain non- orthogonal decomposi- 
tion. We leave this point out, since we concentrate on the analytical aspects of our 
problem. 

We wish to establish the boundary conditions for tp on 9Dq. We will employ the 
boundary condition for u, namely that (u(x),n(x)) = for x e 9D, where n(a;) 
is a normal to 91D) at x. Observe that while in the two-dimensional case, the fact 
that the velocity field u is orthogonal to the normal to the boundary automatically 
implies that it is proportional to the tangent to the boundary, it is no longer true 
in the three-dimensional case. But wc will exploit the helical symmetry and the 
orthogonality condition to compute the tangent vector to 9Do. 

Lemma 2.15. Let D 6e a helical simply connected domain with twice differentiate 
boundary and lef^Q = {{x,yy\{x,y,Qy G P}. Let (a;,t/)^ e 9Do andt = {txity)^ 
be a tangent vector to 9Do (with respect to R^j then t at the point {x,yy is pro- 
portional to 

V = {uxiK^ + y'^) - xyUy, Uy{K^ + x^) - xyuxf , 
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provided that the vector v is not zero. 

Corrolary 2.16 (Boundary conditions for tp). Let Dq G K'^ be a simply connected 
domain with twice- differentiable boundary, then one can choose tjj such that 

(2.60) VlaDo = 0. 

Proof of Corrolary 2.16. By Lemma 2.15, we have on the boundary of Do 



(2.61) 



dt dx ^ dy ^ 

= c(-[(k^ + y'^)ux - xyUy][ux{K^ + y^) - xyuy] 

+ [ux{k^ + y"^) - xyuy\[uy{n^ + x^) - xyux]) = 0, 



where wc used Equations 2.57, which defined The constant c is the proportional- 
ity constant, which depends on the point [x, y)^ . If the conditions of the lemma do 
not hold, namely {ux{k^ +2/^) ^ xyuy^ Uy{n^ + x^) — xyUx) = then Vtp{x, y) = 0. 
Therefore, tp is constant on the boundary of a simply connected domain at every 
point {x,y)'^ G dDo- And since tp is defined only up to a constant, we choose 
V'laDo =0. □ 

Proof of Lemma 2.15. Let n = {hx, ny,nz) be a normal to dH at the point {x, y, 0)^ G 
9Do X {0} C 9D. Since ED is a helical domain, then all the points of the form 
Sp{{x, y, 0)^) belong to 911). Therefore, the field of tangents to the helices ^ is tan- 
gent to the boundary of P. By definition, n is normal to all the tangents to 9D. 
Thus we conclude that 

(2.62) (l',n)=0, 
therefore 

(2.63) yn^ — xuy + kUz = 0. 
This means that n has a form 

(2.64) n= {nx,ny,--nx + -ny)'^. 

Also, by the orthogonality assumtpion, the velocity field u is also orthogonal to the 
helices and thus u also has the form 

(2.65) U= {Ux,Uy,--Ua; + -Uy)'^. 

K K 

But, since u obey the no-flow boundary condition, then (u, n) = 0. We combine 
equations (2.64) and (2.65) with the boundary condition to get 

. y X y X . 

= UxTix + UyHy + (--'"k + ""f/v "'^^ + ~%) 

= i-^ -y) [ ,2 ) J = {u. uy) G [^^^ 



where G is the matrix 

(2.67) G = 



+ y^ —xy 
—xy + 0? 
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Therefore, we conclude that the vector ) is proportional to G ^ ( ]■ Com- 
puting: 



G- 



+ y'^ -xy \ ^ f-Uy 
—xy + x'^J \ Ux 

1 I ^y^ —xy 



(2.68) 



k;2 + a;2 + y2 \^ -xy + x'^ 

1 / 2 2 

i / — K — X —xy 



2/2 \^ -xy -k"^ - y'^J \^ 

/ {k^ + x^)uy - xyux 

/■..2 , ..2\.. 



Observe that wc stipulated that the vector ^ ^"^^ zero in the conditions 

of the lemma. Therefore, the vector {nx,ny)^ is proportional to the vector ((k^ + 
x'^)uy — xyux, xyuy — {k? + y^)ux)'^- Next, let t = {tx, ty) be vector tangent to 9Do 

at {x,y)'^ then the vector {tx,ty,0)'^ is tangent to D at the point {x,y,Q)^. This 
means that {tx,ty,0)"'" is orthogonal to the normal n. From this we conclude 

(2.69) {itx,ty,Of,n)=txnx+tyny = 0. 
Therefore (tx,ty) is proportional to (~ny,nx) which means that 

(2.70) {tx,ty) = c((k^ + y'^)ux - xyUy, (k^ + x'^)uy ~ xyux), 

for some constant c, which depends on {x,y)'^. This concludes the proof. □ 

We now wish to express the connection of the vorticity and the stream function 
and also the dynamics of the vorticity in the new formulation. By Lemma 2.11, we 
can concentrate on the dynamics of cij = fl^. By the same lemma, it is a helical 
function, which we can restrict to a function on TDq, which we denote also by lo. 
Using equation (2.59), we represent ui in terms of tjj: 

dux duy 
dy dx 

1 / -xy -.A /M' 



.>--.)(*).. 

Let us introduce a notation for the matrix-valued function, which appears in the 
previous equation 



(2.72) K{x,y) 



1^ + y^ —xy 
«;2 -I- a;2 -I- y2 y -xy -|- x' 
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We will denote the connection between ui and ^ by an operator Lh- 

(2.73) uj = Ch^- 

With the notation for K, the operator £h can be written: 

(2.74) ChiP = V • {KVxIj) = div(i4:gradV'). 

In the last step, we will write the equation for lj using the stream function formu- 
lation. 

Lemma 2.17 (Dynamics of the vorticity in the stream function formulation). The 

vorticity function u> obeys the following equation: 

(2 75) _ di)_du _^ d^pdu^ ^ 

dt dy dx dx dy ' 

where f^^- 

OX oy 

Proof. Our starting point is Corrolary 2.13: 

(2.76, ^ = 

We have: 

Du) dto dco dco dto 

Olj dui dco 1 duj du 

(2.77) = 97 + + + ^(-^"^ + 9^ + ^9^) 

We now apply the stream function formulation (2.59): 

(2.78) — = — + 
^ ^ Dt dt 



+ - 



1 



,2 



(a± d±\ xy K' + y\ I k- + y -xy \ / ^ 



- dt v^^ ^y) \-i o) \ 



\'By J d^ d^ dy 

Inserting the last equality into equation (2.76) proves the Lemma. □ 

3. Functional Setting and Weak formulation 

3.1. Function Spaces. In this subsection, we will consider measurable functions 
and the appropriate functional space. We will understand helical measurable func- 
tions and helical vector fields as in Definitions 2.2 and 2.4 with the equality being 
true for almost every x. 

We will again use the notation Dq = {{x,y)\{x,y,0) G B}. Observe that the 
assumption on P imply the following for Hq: 

• Bq is simply connected, 

• Do is bounded, 

• Do has a twice differentiable boundary. 
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For a helical integrable function / : 
defined by 



i, the appropriate reduction is / : ©o 



2tt 



(3.1) 



f{x,y) = ^ J f{Sp{x,y,0))dp, 



clearly for g continuous, we will have g = ^IdoxO- The normalization constant 
comes from the computation with / = 1. In what follows, wc will identify / with 
/ and think of / as being a function on Dq. It can be proven that the spaces of 
functions, which we define for function on Dq are equivalent to the spaces of helical 
functions on D but wc will not need this point in our analysis. The only fact, which 
is obvious in our reduction is 



(3.2) 



ess sup / ~ ess sup /, 



which says that the class of csscntialy bounded helical functions on D is equivalent 

to the class of csscntialy bounded on Do. 

Definition 3.1. Denote by LP(I])o) the space of (real-valued) measurable functions 
on Do with the norm 



(3.3) 



\f{x,y)fdxdy. 



Do 



For the inner product in i^(Do) we will use the notation: 
(3.4) if, 9)= j f{x,y)g{x,y)dxdy. 



Do 



Definition 3.2. Denote by W"''^(Do) the space of measurable functions on Dq 
whose distributional derivatives up to order n are Lp(Do) functions with the norm 

d"f ^ 



(3.5) 



IIm'i,p(Do) 



= E 

a|<n 

Also, the space Hq{J])q) is the space of W^'^{I]>o) function, which vanish at the 
boundary, with the norm 



(3.6) 



\Hi = J 

Do 



df 



dx 



+ 



df 



dy 



dxdy. 



We will also introduce the following twisted inner product 
Definition 3.3 ("Helical" inner product). 

{{f,9))h=j ^2j^j.2_^_y2 yox oyjy _xy K^+X 

(3.7) "° 



2 ] ( |2 J dxdy 



= j f{x,y),K{x,y)yg{x,y))dxdy, 



where function K was defined in Equation 2. 72. The appropriate norm is denoted 

by ll/IU = ((/,/))r. 
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Claim 3.4. The norm || • \\h is equivalent to the following semi-norm: 

Proof. The matrix valued function K is symmetric and has eigenvalues 1 and 

2 

-3-7^53—7. Therefore, since Dq is bounded, we have 

K -\-x ~\~y 

(3.9) ;?irR2(^-^'^-^) - (^^Z'^-^) ^ (W,v/), 

where is the maximal distance from in Pq- Integrating over Pq and taking 
square root proves the claim. □ 

3.2. Elliptic Regularity. In this subsection, wc will quote the elliptic regularity 
theory, which was developed by Yudovich in a sequence of papers [21], [22], based 
on the Calderon-Zygmund theory for singular integral operators [18] and Agmon et 
al. [1] results on the regularity of elliptic operators up to the boundary. Yudovich 
[23] applied his elliptic regularity theory to the two-dimensional Euler equations. 
This theory is the main ingredient of our work. We quote his theorem [23, Theorem 
2.1] 

In what follows we will sometimes use .ti,.X2 for the variables x,y. 

Theorem 3.5 (Yudovich's Elliptic Regularity). Let L be an elliptic operator with 
Holder continuous coefficients. Let D be a compact domain in K" with twice dif- 
ferentiable boundary such that L satisfies the following inequality 

(3.10) j u{Lu)dx>m\\u\\]ji(^j^y VueH^{D), m>0. 

D 

Then given f G Lp{D), there exists q G W^'P{D) which solves the problem: 

(3.11) Lq = f, q\an = 0. 
The function q obeys the estimate: 

(3.12) \\q\\w2-P(D) < Cp\\f\\LP{D)- 

Definition 3.6 (Helical operator Ch)- The operator Ch is defined by: 

2 



. . , dx, dx. 



where 

(3.14) K{xx,X'2) = 



+ x\ —X1X2 
2 I —X1X2 + x\ 



Claim 3.7. The operator Ch obeys the following identity, whenever both sides are 
defined 

(3.15) {{f,g))^ = -{CHf,g), 

with f, g vanishing at the boundary. 
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(3.18) J^Hq = ^, g|aDo=0. 



Corrolary 3.8. Let Ch be an operator from Definition 3.6 and f G LP(Do) then 
there exists ip G Vl^^'^'(Do) that solves the following problem: 

(3.16) CHi^ = f ^|aDo = 0, 
and there exists a constant C independent of p such that 

(3.17) llV'lk^.. <Cp||/||p. 

Proof. Clearly, Cr has bounded, smooth coefficients. Combine identity (3.15) with 
Claim 3.4 to verify the inequality in equation (3.10). □ 

We will also need the following statement: (See [23, Theorem 2.2]) 

Theorem 3.9. Let q be a solution of 

Then q obeys the following estimate 

(3.19) \\q\\w^-v<Cp\\f\\L,. 

3.3. The weak formulation and statement of the main theorem. We will 
now formulate the weak problem. As a basis, we use Equation 2.75. Written in 
terms of the stream function tp and operator Ch the equation reads 

We define the weak problem in the following way. 

Definition 3.10. Given F : Do x [0, T) — > and ipo -^o ^ R integrable functions. 
Then ?/> : ©o x [0, T) ^ R with integrable ^ and jCni^ is a solution of the weak 
Euler equation if it obeys the identity 
(3.21) 

T T 

jCHipo{x)(l){x,0)dxdy - J J jCHip{x,t) — {x,t)dxdydt + J j —CHi}-^dxdydt 

Do Do Do 

. } [ fdFy dFA , , , ,^ 

—LHip-^dxdydt = / / - (j)dxdydt, 

ox dy J J \ ox Oy J 

Do Do 

for every smooth function (j) with compact support in Pq x [0)^)- 

We are now ready to state the main theorem. 

Theorem 3.11 (Main theorem). Let F : x [0,T) -> g^^j^ ^/jq^ F,^- 

^ G L°°(Do X [0,r)). Let ipo G W^'P such that £i/Vo e i^(Do) and iZ-olaDo = 
then there exists a unique ip : x [0,T) —^ R with Cnip G -^^(Bo x [0,T)) 
and ipldOo = 0,Vt G [0,T), which solves the weak Euler equation in the sense of 
Definition 3.10. 

Remark. We will also prove that the problem is well posed in Hq norm. 

Remark. Using the ideas of Kato [14], we will prove prove that the solution is 
unique within the class of all stream functions ip with Ch'4> € L°°{[0, T) , i^(Do)). 

We wil prove uniqueness in Section 4 and existence in Section 5. 



18 



HELICAL EULER 



4. Proof of Uniqueness 

In this section, we will prove the uniqueness part of the Theorem 3.11. We will 
employ an energy-principle type of argument, which is identical to that of Yudovich 
[23, section 3]. We follow closely his idea, retaining even some of his notations. 

We will need two lemmas. 

Lemma 4.1. A weak solution of the Euler Equation is a solution of the problem 
for any cylinder Dq x [ii; ^2]; for <ti <t2<T. 

Proof. We substitute into the integral relation a test function of the form h{t)(p{x, t). 

Wc have the following equality 

(4.1) 

T T 

-/i(0) j {jCHil'o){x,y)(t>{x,y,0)dxdy + j h'{t)Mt{il>,^)dt + j h{t)Nt{il>,^)dt = 0, 

Do 

where 

(4.2) MtiiP,^) = - / {CHi^Wxdy, 



Do 



(4.3) iV,(V>,<^) = -j (£h^) ^+CH^{jz'ir. - ¥:'^)dxdy. 



' dx dy dy dx ' 



Mt and Nt are essentially bounded measurable functions. If we choose /i(0) = 
then 

(4.4) j^Mt = Nt, 

in the sense of weak derivatives. Therefore Mt is weakly differentiable and as such 
it is absolutely continuous. Thus we have 

t2 

(4.5) Mt, -Mt,= J NrdT. 

tl 

Substituting the definition of M,N, equations (4. 2), (4. 3) into the last equation 
gives 

(4.6) j jCHtl'{x,y,ti)^{x,y,ti)dxdy - j CHtl'{x,y,t2)^{x,y,t2)dxdy+ 

Do Do 
tl Do 

which means that tp is the weak solution in the cylinder Dq x [ti , ^2] • D 
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Lemma 4.2. Let tl){x,y]t) he a weak solution. Then weak derivatives 
exist and for any p> Po > ^ 



(4.7) max 

0<t<T 



dxdt 



, max 

LP(Po) 0<t<T 



dydt 



< 



LP(Po) 



Proof. ^ satisfies formally the following equations: 



(4.8) 



(4.9) 



+ 



d 



dy J ' dy 



dtp 
dx 



dtp 
'dt 



0. 



Take r > 2 and build a sequence {gk{x,y^t)},k = 1,2... such that Cngu are con- 
tinuously differentiable in Dq x [0,T) , gk\dBo = and CnQk — > in L^"^ = 
L''([0,T];LP(]D)o)) for every p > 0. Define qk by: 



(4.10) 



d_ 

dy 



Cngk 



dgk 
dx 



d_ 

dx 



Cngk 



dgk 
dy 



dFy dF, 
dx 



The application of Theorem (3.9) gives the following bound: 



(4.11) 

Define 

(4.12) 

Then 
(4.13) 



r 9gk 



LP 



+ 



r dgk 

i-Hgk-^ 

dy 



dy ' 



LP 



i'k=i'o{x) + j qk{x,T)dT. 





d^^k ir 

dxdt II LP-'- 



< Clp' milp-r + WjCni^Wloo II . 



The same estimate holds for 

By definition, the solutions of Equation 4.10 obey: 

(4.14) 



- {{lkA))h = j {^Hqk)<pdxdy 

Do 



for any smooth test function (p{x, y, T) 



Do 

0. We transform the relation: 



(4.15) 



-{{Qk,(p))h = 



_//dtPk 
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Wc integrate this from to T to get 
(4.16) j {K{x,y)V^o,V(t){x,y,0)dxdy+ 



Do 



7' 

r I ^'f^ ^ r fdgkd4> dgkd(j)\ [ 1 1 a.\ 

-CH^k^^ + Ch9, [-q^^- ) dxdydt = J (/, cf>) dt. 



Do 

But according to the construction of gk we have 
(4.17) hm / Cngk -r^^ ^ dxdydt = 



k-yoo J \ dx dy dy dx 

Do X [0,T) 



/ 



f dipd(f) dipd(l)\ , , , , 

^hV -^--5 "H-^ dxdydt. 

\ ox ay ay ox J 



Do X [0,T) 

Therefore, by comparing Equation (4.16) with the weak formulation - Equation 
(3.21) and taking the limit, we conclude that: 

T T 







Since ^ are dense in i/^([0, T], ilg) and the sequence f/'fc is bounded in i^([0, T], i?o)- 

Also V'fe arc compact in L'^iPn x [0,T)) because is bounded in W^'f, by Equa- 
tion 4.13 (see [7, Lemma 8.4]). We conclude that ipk converge weakly to tp in 
L'^{[0,T],Hl) and therefore they converge in norm in the space ^^(Bo x [0, T"))- 
Therefore, the derivatives of ijjk converge to limits, which are weak derivatives of 
ij). These derivatives obey the estimate 

(4.19) llSlI^-"- ^ Cp{\\F\\L.,r + \\CHnLA'^nLv,r) , 

where we take the limit r — > oo to establish the required estimate. □ 

Lemma 4.3. Let tpi be a solution of the Euler equation with Cn'tpi € Z/°°(Do X 
[0,T)) with initial condition ipoi with Cni^m G L°°(Do)- Then for every rj > 
there exists 6 > such that if tp2 is a solution of the Euler equation with Ch'4'2 G 
L°°(Do X [0,r)) with initial condition '002 with Ch'4'02 G i°°(Do) then if \\tpoi — 
V'02||h1(Do) < ^ ^^^''^ IIV'i(^) ~ V'2(I')||iji(Do) < V- I'fi' particular, the solutions are 
unique. 

Proof. Form 

(4.20) a = Vi - ^2- 

We use Lemma 4.1 for and t and substract the weak formulations to get 

'^■^i d(j) dtpi d(j)^ 



(4.21) ^ 



dx dy dy dx 



+ £hV'2 - dxdydt = 0, 

\dx dy dy dx ' 
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for any function <p(x, y, t) which is smooth in Dq x [0, T) and (j) =0. 

We now put (f) = a. We have proved that ipi,tp2 and a possess at least 2 weak 
derivatives wliich arc essentially bounded. Therefore, the left hand side of Equation 
(4.21) as a functional on cj) is PF^'^'-continuous and therefore we can substitute 
(f) = a. We get 



(4.22) -{{a,a))^ + 

Do 

We can integrate by parts: 

t 

l\\a\\l= J jK{x,y) 

(4.23) 



at 



dtpi da dtpi da 
dx dy dy dx 



dxdydt = 0. 



Do 



dxdy 



daV 



da 
dy 



+ j j K{x,y) 



Do 



d'^ipi d'^^i\ da da 
dy'^ dx'^ J dx dy 



dxdydt 



dxdydt. 



Since ^ are in W^'^ and thus are bounded, then we can write 



(4.24) 

Denote 

(4.25) 

Differentiating in t, we get: 

dz 

(4.26) < 



\Vaf < M^. 



z^{t) = \\a\\ 



Do 



d^^Pi 



dxdy 



dy'^ 



dx'^ 



Using Equation (4.24): 

(4.27) < CM' 

^ ^ dt - 



/( 



aVi 



dxdy 



+ 



aVi 5Vi 



dy'^ 



dx"^ 



\ya\'^dxdy. 



{Vay-'dxdy. 



We now apply Holder inequality and Yudovich's elliptic estimate to get: 



(4.28) z'^ < CM' 
^ ' dt - 



dxdy 



+ 



d'^tpi d^'ipi 



dy'^ 



dx'^ 



z^-" < M^C-Miz^ 
e 



where Mi = ||-ChVi||l°°(Dox[o,t))- We integrate 

(4.29) z\t) - z\Q) = {CMy2CMit. 



(4.30) 

Choose ry > 0. Denote 
(4.31) 



CM 



< 2CMit ■ 



£(0) 
CM 



K 



ACMi 
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We will chose a sequence 61,62, ■■■Sk+i in the following manner. Take 6k+i = rj 
and with 5i+i, ..,6k+i chosen, choose Cj such that 



(4.32) 



With this choice of Ci, choose 6i such that 



CM 



< 



1 



41/ei • 
T 



(4.33) 
Denote 

^'■''^ ■ ACM, 

For this sequence of choices, the following holds: if z{[k— l]r) < 5% then z (kr) < 
di+i. This is because using e, we apply equation (4.30) starting from time {k — 1)t 
then 



z{kT) < (2CM1T + 



(4.35) 



z{[k - 1]t) 



CM 



< 



<6. 



i+i- 



Therefore, if wc have two initial conditions such that ||'0oi ~ "00211//! < i5i then 
||0i(T)-V'2(T)||hi Inparticularz(i) =0 for z(0) = 0, for every te [0,T). □ 

We will now prove uniqueness in a larger class of functions, following the work 
of Kato [14]. Wc will need the following Sobolcv-typc inequality (see [12, II. 2. 2]). 

Lemma 4.4. Let (j) in Hq{D) then for 2 < p < 00 

p-2 

(4.36) UWl. < ' ^ 



V2V2 



Proposition 4.5. Let tpi,tp2 be a weak solution of the Euler Equations such that 

Cui'i e L°°(Do X [0,r)) and Ch'4'2 e L°°([0,r),L2(Do)) that have the same 
initial condition tpo with £//V'o G L°°{V)o) , forcing F e L°°{3o x [0,T]) and 



/ = ^ - ^ e L-(Do X [0,T]) then 

(4.37) ^1 = V'2. 

Proof. Continuing along the lines of the previous proof, form the difference 

(4.38) a = Vi-02- 

Again, use Lemma 4.1 for and t and substract the weak formulations to get 

(4.39) { •' ^ 
+ -C//V2 ( - 



dx dy dy dx 
dxdydt = 0, 



/ da d(p da d<j) 
\dx dy dy dx 

for any function (t>{x, y, t) which is smooth in Dq x [0, T) and 



We now put 



0. 



a. We have proved that tpi,ip2 and a possess at least 2 weak 
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derivatives that are in L^. Therefore, the left hand side of Equation (4.21) as a 
functional on (j) is ff^-continuous and therefore we can substitute (j) = a. We get 



at 



(4.40) -((a,a))^ + 
We can integrate by parts: 



da dipi da 
dx dy dy dx 



dxdydt = 0. 



(4.41) 



Differentiating in t 
1 d 



Do 
t 



da 
dx 



da 
dy 



K{x,y) 



Do 



dxdy 

d'^ipi d'^ipi \ da da 
dy"^ dx"^ J dx dy 



(4.42) 



K{x,y) 

Do 

+ j K{x,y) 



da 
dx 



da 
dy 



dxdy 

dSh d'^ipi\dada 
dy"^ 



dx'^ J dx dy 



dxdydt 
dxdydt. 

dxdy 
dxdy. 



Since Cni^i € L°° then we can apply elliptic estimates 



(4.43) 



dxidxi 



< CpM, 



LP 



where M = sup ||£ij'^i||j,oo(]D)). We use Holder inequallity to get 

te[o,T) 



(4.44) 



^^||af,<CpM||Va||i.„ 



where i + i = 1. Since the helical norm \\a\\h is equivalent to ifg-norm ||Va||i2 
and since ajoDo = we can rewrite the previous incquallity as 



(4.45) 



ld_ 
2di 



\a\\l < CpM\\Va\\l., 



We use Lemma 4.4 to conclude 



(4.46) 



llVal 



L29 



< 



{V2 



2(i-i/g) 



l|Vaf//||Val|^^^/^\ 



Observe that by applying elliptic estimates for p = 2, we have 

(4.47) llVallni < 2C\\Ch^ji - ChH\lhb„), 

which is bounded in time by the assumptions of the proposition. Denote 

(4.48) z = \\Va\\h. 
Collecting the estimates and redefining C, we have 



(4.49) 



V2) 



2(1-1/9) 



.2/9 _ 



= Cp 



p 



{p-l)V2 



2/p 



^2-2/p 
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(4.50) ^^z'/P <Cp( 

^ ^ 2dt - ■^V(p-l)\/2 

Therefore, 

(4.51) z{t) < P M ctyi''. 

[p- l)\/2 

Taking t < and jj — > oo proves that z{t) =0 and by iteration, for all times. □ 

5. Proof of Existence 

In this section, we will prove existence of the weak solutions. Our strategy is to 

smoothen the initial data and the forcing, apply classical methods to prove exis- 
tence of the solutions for smoothened data and then apply compactness arguments 
to extract a subsequence converging to the weak solution. A crucial step in the 
argument is an application of Lemma 4.2, which allows the required control of the 

compactness. 

5.1. Construction of a sequence of smooth approximations. 

Claim 5.1. Let C € R" 6e a domain with twice differentiable boundary, with 
compact closure, let f G L°°{C). Then there exists a sequence fn G C°°{C) such 
that fn{x) converges to f{x) for almost every x G C and ||/n||L~ < 

Proof. Let a -.W^ —^ [0, oo) be a standard mollification kernel (A positive, C°°(M") 
function, which is compactly supported on the unit ball and whose total integral is 
one) . And define 

(5.1) -e(.) = Xf). 
Define the set by 

(5.2) Ce = {x€ C\d{x, dC) > e], 

the set of all a: € C whose distance from the boundary of C is greater then e. Take 
a sequence of £„ going to zero and define 

(5-3) fn = (Te„ * {xC2.„f ), 

where * is the convolution in M". The sequence {/„} obeys the requirements of the 
lemma. □ 

Observe that by employing Lebesgue dominated convergence, we conclude that 
fn converges to / in every LP for p < oo and weak-* in L°°. 

Let us now return to the weak problem as defined in Definition 3.10. Let oJo = 
Ch^Po be the initial vorticity and / = — be the forcing of the vorticity. 
We choose the sequences e C°°(Do) and /" e C°°(Do x [0, T)), which obey the 
conditions of the Claim 5.1. Define JIq = -^^t^o ^^'^ k^-^"' ^^'^ vector field ^ 

was defined in equation (2.6). We know extend fig and f" to helical vector fields 
on D by Equation 2.51. 

We now obtain a sequence smooth initial velocity fields Uq and body forcing 
terms F", which obey the orthogonality condition. We now wish to solve the Euler 
Equations in the vorticity formulation with initial data Uq and vorticity forcing 
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term F" in the domain D . We will employ the following theorem of Ferrari^ (see 
[11, Theorem 2]). 

Theorem 5.2 (Ferrari). Let Hq G H^{D),s > 2 be the initial vorticity , D is 
bounded simply connected domain then there exists a solution ft G C{[0,T],H^{D)) 
if and only if 

T 

(5.4) J \\u)\\Loo(^o)dt is finite. 



Theorem 5.2 is an extension of the theorem by Beal, Kato and Majda [4], who 
proved the same criterion for Euler equations in the whole M". We will need a 
modification of Theorem 5.2, for a helical domain ED, which is bounded in the x, y 
direction and periodic in the z direction and solution, which are periodic in the z 
direction. It is proved along the same lines, as the original theorem of Ferrari. We 
now apply the theorem of Ferrari. 

Lemma 5.3. For every n S N there exists a smooth helical solution of the Eu- 
ler Equations fi" € C°°(D x [0, T)) with the initial vorticity CIq and forcing in- 
Moreover, it obeys the following estimate: 

(5.5) sup ||n"(t)||ioo < M, 

te[o,T) 

where 

T 

(5.6) M = C{\\no\\Loo + J \\f{t)\\Loo). 



Proof. Observe that if Jl" obeys the required estimate then it fulfills the conditions 
of the theorem of Ferrari and it is in the class H". Then we take s to infinity 
will prove that 17 is a C°° vector field. Observe that ft are classical diffcrcntiable 
functions and therefore, we can apply the standard tools of calculus. We will use 
Corollary 2.13 to conclude the estimate for fl^. Then we will use Lemma 2.11 and 
the fact that B is contained in a cylinder of finite radius and therefore \Ctx\ and 
\Qy\ are smaller then C\flz\ for some bounded constant C. □ 

From now on we will fix p > 2. Unlike in the previous section, where the 
dependence of the constants on p was crucial to the argument, in this section it 
will not play any role, since p is fixed and we will subsume this dependence into 
constants. 

After obtaining the solutions r2"(x, t), wc define 

(5.7) a;" = (fi"),, 

a sequence of smooth helical functions, which we reduce to function on Dq- Then 
we obtain solutions to the problem 

(5.8) CHr=co", r\om = (^- 



■We restate the theorem adapted to our setting in a standard manner. 
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5.2. Uniform bounds on approximating sequence. 

Claim 5.4. Let ^" he the stream function corresponding to then tp" obeys the 
following estimate: 

(5.9) sup||V"W||v^2.p(D„)<CM. 

Proof. This is the immediate corollary of the elliptic regularity, Corrolary 3.8 to- 
gether with the fact that the L°° norm of the vorticity is uniformly bounded. □ 

We will also need a uniform bound on the time derivative of tjj"': 
Lemma 5.5. Let tjj" be the stream function corresponding to 

(5.10) sup < C{\\fU^ + M'). 
te[o,T) 

Proof. We will apply Lemma 4.2. The estimate in Lemma 4.2 depends on the 

norm of the body forcing function F". Wc wish to replace it with the estimate 
on f". Since the velocity fields are incompressible, we can replace F" with it's 
incompressible part without affecting the dynamics. But then F" obeys 

(5.11) AF" = VAf". 

We now apply the elliptic regularity in the form of Theorem 3.9 to conclude that the 
norms of the body forcing term are controlled by the vorticity forcing term. Then 
we can apply Lemma 4.2 and the uniform estimate on ||V'"(*)||w2,p from Claim 5.4 
to conclude the lemma. □ 



5.3. Extraction of a solution to the weak problem. After establishing uniform 

in n estimates, wc will now apply Aubin compactness argument for the sequence 
{■0"}. We will use the following theorem, [7, Lemma 8.4] 

Theorem 5.6. Let Xi,Xo,X-i be reflective Banach spaces such that Xi is com- 
pactly embedded into Xq and Xq is continuous embedded in Let Um be a 
bounded sequence in LP^{[0,T),Xi) such that ^ is bounded in LP^{[0,T) ,X_i) 
then there exists subsequence Um' which converges in L^i([0,T) ,Xo). 

Corroleiry 5.7. There exist a subsequence ofip"' which converges in 

Lf([0,T),W^i'f(Do)). 

Proof We will apply Theorem 5.6 with Xi = W^'P(Po) and Xq = X^i = W^-p{I3o). 
The uniform estimates are Claim 5.4 and Lemma 5.5. Observe that we obtained 
estimates in L°°([0, T) , Xi) for the functions and L~([0, T) , X_i) for the time 
derivatives, clearly they imply LP([0, T) , ATi) and Z/^'([0, T) , X_i) estimates, re- 
spectively. □ 

Proposition 5.8. Let tpo be such that Cni^o = ojo € L°°(I])o) and F G L°°(Do x 
[0,T)) then there exists tp such that Cni^ = w e Z/°°(Do x [0,7")) which obeys the 
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identity 
(5.12) 



u)o{x)(j){x,0)dxdy ~ J j Lo{x,t) — {x,t)dxdydt + J J -g-co—dxdydt 

Do Do Do 

dtp d(j) , , , } f fdFy dF^ 



—uj—dxdydt = [ - ^— (pdxdydt, 
dx dy J J \ dx dy J 

Do Do 

for every test function (j). 

Proof. Let Cli" = O" bo the sequence of functions, which wc built in Subsection 
5.1. By Lemma 5.3, they are uniformly bounded in L°°(Do) and therefore, we can 
extract a subsequence, which converges in weak-* topology to a function u). By 
extracting from the subsequence a subsubsequence, we see that will converge 
in L'''([0, r) , M^^'^(Do)) to a function ijj. By uniqueness of the limits, V' must be 
the stream function which corresponds to uj. Therefore VV'" will converge to V'tp 
in i:P([0,T),LP(Do)) C LP(Do x [0,T)). Therefore cj"^ converge weakly in 
LP(ID)o X [0, T)). For any smooth function;/) with compact support in Dq x [0, T),the 
functions ipn is the weak solution of the Euler Equations in the sense of Definition 
3.10 and therefore they obey 

(5.13) 

T T 

j ijj'${x)4>{x,Q)dxdy - j J ij"^{x,t)dxdydt + J J ^^u^'^dxdydt 

Do Do Do 

T T 



1 1 5^."||^*.< .jl(?S..?S.^ 



Do Do 

We now see that all the terms in this equation converge weakly to the appropriate 
limits and therefore, t/j, the limit, will obey the required equality and therefore it 
is the solution of the Euler Equations. □ 
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